1 Introduction

Coastal Resilience in West Africa and the nearby Islands of Cabo Verde is challenging because of mounting human drivers associated with population growth, socio-economic development, natural hazards, anthropogenic climate change, and limited observing systems. An estimated 100 million persons live along coastal West Africa and expected to increase throughout the twenty-first century rapidly. Most of the region has a limited observing network along coastlines and also over the oceans. Ship reports provide most of the real-time ocean and weather conditions between Cape Verde and Coastal West Africa, while the Prediction and Research Moored Array in the Tropical Atlantic (PIRATA) of buoys are sufficiently far from coastal zones in West Africa [1]. There are limited tide-gauges and coastal weather stations along coastal West Africa. Automated weather stations in Cape Verde were critical to observing coastal flooding from Hurricane Fred, as well as from sea-level rise and natural hazards [2].

During the northern hemisphere summer season, African Easterly Waves (AEWs) move across West Africa every 3–5 days and may be associated with flooding inland but can evolve further to tropical disturbances over the Atlantic Ocean [3, 4]. Sometimes, an area of low pressure with a well-defined area of vorticity between 850–700 hPa can develop after the AEWs exit the coast leading to a tropical depression or a named storm. They can then be the cause of loss of human life and material when they deepen near the coast. In 1999, the disturbance associated with Tropical cyclone Cindy was observed from Senegal and lead to more than 100 fatalities along the coast of Senegal and Mauritania [5]. Senegal is located in the Western Sahel, with a growing population of approximately 15 million residents. However, many of the residents live near the capital of Dakar, which is the westernmost city along the African Coast (14.59 N, 17.5 W). More than 20% of the Senegalese population lives in Dakar, although its area does not cover more than 0.5% of the country, according to the latest report from the National Agency for Statistics and Demography of Senegal (http://www.ansd.sn/). Senegal has a coastline of 706.72 km long [6] and has located in the southern part of the Canary Upwelling system (a typical Eastern Boundary Upwelling System (EBUS). The Senegalese upwelling season, which generally starts in October and ends in late May to early June, ensures a high level of primary production [7,8,9,10,11]. The Southern Senegal marine ecosystem is well known to be a spawning and nursery area for small pelagic fish species. Therefore, artisanal fishing is the first economic activity generating income in Senegal.

Two types of swells characterize the Senegalese coast: (1) In the northern zone (between St-Louis and Dakar), the trade winds and disturbances of the North Atlantic basin generate swells of north/northwest direction throughout the year. (2) In the southern zone (from Dakar to Cap-Skirring), swells are less frequent because northwest Atlantic swells are diffracted over the peninsula of Cape Verde (defined as the Region of Dakar) which considerably weakens them. However, swells occur mostly in NH summer period and can be generated by Mesoscale convective systems and/or African Easterly Waves which can deepen into tropical depressions off the coast [2, 5, 12, 13]. The consequences can be disastrous for the populations living along the littoral. Indeed, this zone is characterized by a broad and shallow continental shelf (i.e., 20–30 m over tens of kilometers) compared to the northern zone [8, 14]. The shallow depth of the ocean and the concave shape of Senegal's coastline can lead to the gradual advance of the sea, thereby increasing coastal population vulnerability to flooding during strong swell events.

Hurricanes rarely occur in the extreme Eastern Atlantic Ocean, but on 31 August 2015, Hurricane Fred passed through the Cape Verde Islands, causing damage in seven of the ten Islands [2]. This storm rapidly intensified from a tropical depression at 0000 UTC on 30 August to a hurricane 0000 UTC 31 August. However, prior to reaching Cape Verde, this storm caused seven fatalities associated with a fishing vessel and caused damage along Senegal's coastline. On 30 August 2015 reports by administrative authorities of Fatick region, which is in southwest Senegal, reported significant losses in the villages (Palmarin, Sessene, Dionewar, Ngallou, Diakhanor, Djifere, Cap Skirring) caused by a strong swell during the night of 30 August. They report that approximately 20 houses were destroyed, displacing many people. Further north, media reports suggest that more than 200 homes were damaged/destroyed along with boats. Residents that were interviewed in Bargny, Senegal, and the Hann District of Senegal spoke of the rise of the seas during the day of 30 August, leading to widespread destruction (Figure S1).

The objectives of this study are to use a set of fine grid spacing simulations using the Weather, Research and Forecasting (WRF) model to characterize the disturbance associated with TC Fred through the period of 29 and 30 August during its evolution. We then use WRF output from these simulations to estimate Ekman transport of water toward the coast of Senegal that could have been associated with coastal damage on 30 August 2015. This paper is organized as follows: the used data, model simulations and Ekman transport estimates are presented in Sect. 2. The results and conclusions are presented in Sects. 3 and 4.

2 Data and methodology

2.1 Data

Based on satellite observations and model analysis, the National Hurricane Center (NHC) maintains a climatology of all Atlantic tropical cyclones since 1851 in the whole Atlantic and the Pacific, called HURDAT. For each storm, HURDAT contains estimates of the latitude, longitude, 1-min maximum sustained surface winds, minimum sea-level pressure with a precision of 0.18° latitude/longitude, 2.57 m/s, and 1 mb, respectively [15] for every synoptic time (0000, 0600, 1200, and 1800 UTC). Here, the six-hourly positions are used to track the evolution of Hurricane Fred. These data were downloaded from the NHC website (https://www.nhc.noaa.gov/). The new gridded Daily Advanced Scatterometer (DASCAT) wind and wind stress data are used to examine 10 m wind associated with Fred genesis. Data are estimated over global oceans from Advanced Scatterometer Data Products (ASCAT) retrievals using objective methods and have a spatial resolution of 0.25° in longitude and latitude. The calculation of daily estimates uses ascending/descending tracks available and valid retrievals. According to the ASCAT sampling scheme, the objective method allowing the determination of regular in space and surface wind fields uses ASCAT observations as well as ECMWF analyses. The latter is considered as the temporal interpolation basis of ASCAT retrievals [16]. More details about data, objective method, computation algorithm may be found in [17].

Satellite Brightness temperature data of the thermal infrared channel (10.8 mm) is produced by the Climate Prediction Center (CPC), NOAA National Centers for Environmental Prediction. The data are provided in netCDF-4 format on a 4-km-equivalent latitude/longitude grid over the latitude band 60°N-S.

ERA‐Interim is a global atmospheric reanalysis produced by the European Centre for Medium‐Range Weather Forecasts (ECMWF). It is produced with a sequential data assimilation scheme, advancing forward in time using 12‐hourly analysis cycles. In each cycle, available observations are combined with prior information from a forecast model to estimate the evolving state of the global atmosphere and its underlying surface. The spatial resolution of the data set is approximately 80 km (T255 spectral) on 60 levels in the vertical from the surface up to 0.1 hPa [18]. The new high-resolution Reynolds sea surface temperature (SST) data is used to analyze the SST pattern associated to Hurricane Fred. The data are stored on a spatial resolution of 0.25-degree and a temporal resolution of 1 day. The information about the interpolation method used to conduct is described in [19].

The winds measured at Dakar Yoff (located at 14° 44′ N, 17° 30′ W, 27 m above ground) are representative of the winds over Southern Senegal [8, 9]. They are also used to evaluate the Weather Research and Forecasting (WRF) model's ability to reproduce the surface wind circulations over the southern Senegalese coast. The tide gauge data of Dakar station from 1993 to 2019 is used to analyze the sea-level rise during the passage of Hurricane along the Senegalese coast. Data are retrieved from http://uhslc.soest.hawaii.edu/data/?rq

The variables presented above are used to analyze the atmospheric conditions occurring during the genesis of Hurricane Fred, but also to analyze the capacity of the WRF model to reproduce the event.

2.2 Methodology

Recent case studies based on observations have shown that Ekman transport can have a strong effect on coastal flooding especially when the storm moves slowly, the storm-induced winds are parallel to the coast [20,21,22]. Ekman transport has been cited as one of the key mechanisms causing subtidal water level fluctuation inside the Chesapeake Bay for a wind forcing longer than 4 days [23]. Based on a 3D storm surge model of the Chesapeake Bay, Shen and Gong [20,21,] studied the impact of model domain size, wind directions and Ekman transport on the storm surge of the extratropical cyclone Ernesto (2006) in the Chesapeake Bay. Their results show that the Ekman setup along the coast contributes significantly to the water-level variations during storm events. It enhances (reduces) surge inside the Bay under the wind forcing from N and NE (SW, S, and SE) directions. By using a set of data from an ADCP station placed at a depth of 65 m on the sea bed off the coast at Sète, supplemented by the ocean circulation model SYMPHONIE, Michaud et al. [22] studied the hydrodynamics in the Gulf of Aigues-Mortes (NW of the Gulf of Lion) during the storm of 18 February 2007. They showed that “the duration of the storm (about 36 h at the apex) explains the continuous increase of the current as predicted by the theory. The frictionally induced Ekman transport explains the current shear in the surface layer in the rising stage of the storm, and the addition of high waves and strong wind at the apex is more in favor of strong vertical mixing in the surface layer.”

To examine the role of Ekman Transport, we use the Weather Research and Forecasting (WRF) model output, which has an outer 12 km grid spacing domain and an inner 4 km grid spacing domain. We use the WRF version 3.71 [24] in this study to examine the tropical disturbance off the coast of West Africa. The National Center for Environmental Prediction (NCEP) final analyses FNL is used to initialize and provide boundary conditions using FNL analyses for the parent domain at 12 km every 6 h. Three simulations using two domains are used in simulating TC Fred. The simulations are started at different times beginning on 28 August 0000 UTC and initialized every 12 h (Table 1). The coarse grid has 12 km grid spacing with 31 vertical levels, 251 east–west grid points, and 181 north–south grid points covering much of West Africa and the extreme Eastern Atlantic (Fig. 1a). The fine grid uses 4 km grid spacing with 31 vertical levels, 196 east–west points, and 160 north–south points covering Senegal and extreme Eastern Atlantic (Table 1). The fine grid simulations begin 6 h after the coarse grid simulation is initiated. There is two-way feedback used between the coarse and fine grid domains. The larger domain uses the Kain-Fritsch cumulus parameterization [25], while the inner domain does not use a cumulus parameterization; Thompson microphysics [26] is used in the outer and inner domains.

Table 1 WRF simulations undertaken for this study, with grid spacing, start and end times
Fig. 1
figure 1

a Map of West Africa including the simulation domains and the orography. The red larger represents the domain 1 and the smaller red rectangle is for domain 2. b Bathymetry (in dashed line) and location of observed damages (in red point) during the passage of the Hurricane Fred along the Senegalese coast

Surface stresses (\({\tau }_{x}\) and \({\tau }_{y}\)) are computed from WRF simulated ten-meter wind components and used to estimate the zonal and meridional Ekman transport with a focus on 30 August 2015 at various locations along the Senegalese coast (Fig. 1b). We computed Ekman transport (Q, m3·s−1·m−1) using a constant atmospheric density (ρa = 1.2 kg m-3), a dimensionless empirical drag coefficient (Cd = 1.3 × 10–3; [10], the empirical formulation of wind stress, a constant sea water density (rho = ρ = 1000 kg m-3) and the Coriolis parameter f as follows:

$${\tau }_{x}={\rho }_{a}\bullet {C}_{d}\bullet \sqrt{{u}^{2}+{v}^{2}}\bullet u$$
(1)
$${\tau }_{y}={\rho }_{a}\bullet {C}_{d}\bullet \sqrt{{u}^{2}+{v}^{2}}\bullet v$$
(2)
$${Q}_{x}=\frac{{\tau }_{y}}{f\bullet \rho }$$
(3)
$${Q}_{y}=\frac{-{\tau }_{x}}{f\bullet \rho }$$
(4)

In these two formulations above, the Qx and Qy correspond to the zonal and meridional Ekman transport components, respectively.

In order to assess the effect of Ekman transport during the period 30–31 of August on the Senegalese coast, we calculated Accumulated Ekman Transport (AET) which is the time-integrated Ekman transport using the following formula

$${\mathrm{AET}}_{x}={\int }_{{T}_{1}}^{{T}_{2}}{Q}_{x}\bullet \mathrm{d}T$$
(5)
$${\mathrm{AET}}_{y}={\int }_{{T}_{1}}^{{T}_{2}}{Q}_{y}\bullet \mathrm{d}T$$
(6)

where \({\mathrm{AET}}_{x}\) and \({AET}_{y}\) correspond to the zonal and meridional Accumulated Ekman Transport components, respectively. T1 and T2 are the initial and the final dates, respectively.

3 Results

3.1 Environmental conditions of the genesis of Fred

On 28 August 2015, a mesoscale convective system embedded in an AEW trough from Mali moved into southeastern Senegal (Fig. 1a, b). It moved slowly, crossing the Guinea Highlands (located around 11° N,12° W), and arrived at the coast between 29th 1200 UTC and 30th 0000 UTC. During this time, there was a closed cyclonic circulation at 700 hPa. Convection gradually strengthened during the day of 30 August, with AEW remaining almost stationary near the Senegal and Guinea coasts. From 30th to 31 August, the convection increases in scale, and there is a strengthening vortex at 700 hPa (Fig. 2d–f). Observations show that the surface cyclonic circulation began to form during the same period, and the system was classified as a tropical depression by the NHC about 0000 UTC on 30 August. The wind sensor of "Melax the buoy" (only Senegalese Buoy) located at 400 km north (14° 20′ N,17° 14′ W) of the disturbance also observed the effects of the tropical depression with surface wind speed records of 12 m/s (about 24 knots), associated with the intensification of the disturbance (figure not shown).

Fig. 2
figure 2

Meteosat Second Generation brightness temperature (in color, units in Degree C), Era Interim 700 hPa wind (in streamline) and positive relative vorticity (contours greater than 1.10–5 s−1) showing the evolution of the disturbance that gave rise to the Hurricane Fred off the Senegal-Guinean coast

The disturbance developed in a favorable environment with ocean surface temperature higher than 28 °C (Fig. 3a) and low vertical wind shear (figure not shown). During 30 and 31 August, the system intensifies further to reach a category 1 hurricane on the Saphir-Simpson scale before reaching the Cape Verde Islands, a situation that was never encountered during the twentieth century with some evidence suggesting a hurricane in 1892 [2].

Fig. 3
figure 3

a Mean daily Sea Surface Temperature for 30 August. Units in Degree C. Data are from https://psl.noaa.gov/data/gridded/data.noaa.oisst.v2.highres.html. Mean daily 10 m wind (color and vector, in knots) in b for 30 August and in c for 31 August

Figure 3b, c shows DASCAT wind fields on 30 August and 31 August. Southerly winds occur near the Senegalese coast on 30 and 31 August. Winds of more than 25 knots (~ 12.5 m s−1) are persistent during this period. The southerly fetch extends for more than 500 km for both days, which could increase the wave heights beyond 4 m on the coastlines to the south of Dakar, especially during the morning hours of 31 August.

3.2 Model evaluation

Tracks of the simulated disturbance relative to the NHC best track are presented in Fig. 4. The black solid line shows the NHC evolution and track of tropical disturbance associated with Fred making the closest approach to West Africa near the Guinea-Bissau at 0000 UTC 30 August. Throughout the day, TC Fred is located several hundred kilometers off the coast of Senegal, and by 0000 31 August, the center of then Hurricane Fred is located approximately at 4.2 degrees west of Dakar, according to NHC archives. In all of the simulations, the forecasted location of the storm is behind the official NHC location. Simulation D2A, which began on 28 August, shows a track that is closest to the NHC track. However, its position at 0000 UTC 31 August is to the southeast of the actual location. It is approximately 6 h behind the NHC position and located near the estimated NHC position at 1800 UTC, 30 August. The track of simulation D2B, which began at 1800 28 August, is displaced to the east of the NHC official track and located much closer to the Guinea-Bissau and Senegalese coastlines. It is also displaced in time and occupies a latitude similar to NHC. The third simulation, D2C, began at 0600, 29 August, and is displaced to the west of the NHC position. It also occupies a biased position and is located at 14° N, similar to the other two simulations at 0000 31 August. The differences in the simulated position of TC Fred, provided a useful means of evaluating coastal winds and the Ekman transport at the Palmerin, Rufisque, and Thiaroye Senegal, which are at or near locations where damage occurred (Figure S1).

Fig. 4
figure 4

Evolution and track of tropical disturbance associated with Fred from 30 August, 2015 at 0000 UTC to 31 August, 2015 at 0000 UTC. Black solid line is NHC data, blue dashed line for simulation 2A, red dashed line for simulation 2B and green dashed line for simulation 2C

Figure 5a–f represents 10 m wind on 30 August 2015, respectively, at 0600 UTC and 2000 UTC. Simulated disturbances at 0600 UTC 30 August are located between 11 and 12.5 N in the three simulations (Fig. 5a, c, e). The simulated center of the disturbance in D2B is located closest to the coast (12.3 N, 17 W), while the center in D2C (11.7 N, 18 W) is located the furthest out to sea; the center of the disturbance in D2A is located further south near 11.3 N. Winds of 30–40 knots are simulated in all simulations but cover a larger area in simulations D2A and D2B. At 2000 UTC, the centers of the simulated storms are located near 13.5 N, 20.5 W (D2A), 14.5 N, 20 W (D2B), and 13.5 N, 21 W (D2C) (Figs. 5b, d, f). The strongest winds of 40–50 knots are found near the center of the simulated disturbances. The three simulations show strong southerly winds (20–30 knots) along the Senegalese coastlines at 2000 UTC. In the regions where damages were reported, southerly winds would drive waves northward toward the Senegalese coast.

Fig. 5
figure 5

Simulated 10 m-winds (color and arrows in kt) on 30 August 2015, respectively, from left to right at 0600 UTC and 2000UTC. Top for simulation 2A, middle for simulation 2B and bottom for simulation 2C

Figure 6a, b shows the observed and simulated meridional wind component at selected locations along coastal Senegal and WRF meridional winds averaged over the nearest grid point of Dakar observation. It should be noted that, except for Dakar, all the other stations present discontinuities synonymous with a lack of data (Fig. 6a). At the southernmost station of Cap Skirring, a rapid transition to southerly winds is found on 30 August, beginning at 0600 UTC. Stations located further north also show a pattern of southerly winds after 0600 UTC at Banjul, Gambia, Dakar, and Ziguinchor Senegal. The wind speeds of 10–15 knots are found at Cap Skirring after 0600 UTC, while winds of 10 to 15 knots are found after 1200 UTC in Dakar, Senegal. This pattern of observed southerly winds between Cap Skirring and Dakar, Senegal (distance 142 nautical miles or 263 km) would have created fetch that could have led to significant wave heights in coastal zones based on the Sverdrup-Monk-Bretschneider Nomogram (Kelvin Rodolfo [27].

Fig. 6
figure 6

a observed Meridional wind (kt) at Dakar (red) and Cap Skiring (yellow) from 29 August, at 00:00 UTC to 31 August 29 at 00:00 UTC, b simulated WRF meridional winds averaged over the nearest grid point of Dakar observation and during the same period as in a blue dashed line for simulation 2A, red dashed line for 2B and green dashed line for 2C. The solid black line represents the observed meridional wind at Dakar

Simulated winds (Fig. 6b) in Dakar also show a pattern of southerly winds in agreement with observations. Simulation D2A shows continuous southerly winds beginning at approximately 0500 UTC, 30 August, and continuing throughout the day. The winds reach the maximum value of ~ 18 knots between 1600 and 2200 UTC. Simulation D2B shows a short period of southerly winds but produces the strongest simulated meridional winds of more than 19 knots between 2100 and 2300 UTC on 30 August. While the simulated disturbance is located furthest west relative to the other simulations, it simulated southerly winds beginning the earliest (0400 UTC on 30 August) and increases throughout the day, reaching its largest values between 1700 and 2200 UTC.

3.3 Zonal Ekman transport estimation

The Ekman transport, which is determined from simulated zonal and meridional wind components at grid-point locations near Palmerin, Rufisque, and Thiaroye (see Fig. 1b), is estimated using the Eqs. (3, 4). Zonal Ekman transport indicates increasing eastward transport toward the coastline after 1200 UTC, 30 August, based on D2A and D2B (Fig. 7a–c). The largest transport occurs after 1800 UTC and is found at Rufisque relative to the other coastal locations in these simulations. Zonal transport estimated by simulation D3C is smaller relative to D2A and D2B. However, all of the transport is eastward in this simulation, with the largest values found around 0000 30 August and steady westward transport between 0900 and 2000 UTC 30 August.

Fig. 7
figure 7

Computed Zonal Ekman transport (m3· s−1·m−1) at grid-point locations near Palmerin, Rufisque, and Thiaroye (see Fig. 1b)

Figure 8 shows varied Accumulated Ekman Transport (AET) between 30 and 31 August from the three simulations. The accumulated Ekman transport computed from simulation D2A occurs to the south of Senegal (Fig. 8a), while a higher accumulated Ekman transport is found in D2B and D2C along the Senegalese coastline (Fig. 8b, c). The largest accumulated Ekman transport along the Senegalese coastline occurs in 2B, which has to closer track to the coast.

Fig. 8
figure 8

Zonal component of the accumulated Ekman transport in color from 30th to 31st august for 2A (top), 2B (middle) and 2C (bottom). Total accumulated Ekman transport is represented in vectors

Figure 9 shows the temporal series of the tide gauge data of Dakar station (longitude = −17.42°W, latitude = 14.67°N) from 1993 to 2019. It reveals the seasonal variability of the sea level due to the southern Senegalese Upwelling variability [9]. However, the maximum height measured during this period is 1.35 m. This value corresponds to the date of passage of Hurricane Fred.

Fig. 9
figure 9

Temporal series of the tide gauge data of Dakar station (longitude = −17.42°W, latitude = 14.67°N) from 1993 to 2019

4 Conclusion

During the evening of 30 August and likely into the early morning hours, residents along the coast of Senegal from Dakar through Palmarin experienced coastal flooding and infrastructure damages associated with the tropical disturbance that became Hurricane Fred. The flooding occurred when the tropical disturbance was developing but several hundred km away from the coast. Discussions with residents of Bargny (located near Rufisque) spoke of rising seas throughout the night and early morning hours, which caused the damage. Additional damage in the surrounding areas of Rufisque may have occurred (Ndoye, personal observation) but not reported. In this work, we used the WRF model to examine the potential impact of the tropical disturbance with a focus on 30 August 2015. The three simulations show a spread in the storm tracks relative to the NHC official track, with simulation D2A, which was initiated on 28 August 0000 UTC having the smallest track errors. Simulation D2B had an eastward bias, while simulation D2C had a westward bias.

Observations between Dakar and southern Senegal and including the Gambia show southerly winds during the morning through the night of 30 August, with model simulations showing agreement with the coastal observations. Satellite overpasses show winds greater than 15 knots along coastal Senegal on 30–31 August, which would likely create a long fetch. The southerly winds may have impacted the coast of Senegal through eastward Ekman transport. High winds, and accumulated Ekman transport produces a storm surge, which is usually well above the high-tide level. The use of a Munk-Sverdrup-Bretschneider nanogram (a diagram used to estimate wave height based on hours of fetch and wind speed) suggests that a fetch of more 400 km for more than 24 h with wind speeds of 20 knots would produce wave heights of more than 2 m along coastal zones south of Dakar, Senegal. The largest accumulated Ekman transport along the Senegalese coastline occurs in 2B (see Fig. 8), which has to closer track to the coast is in agreement with the high value of sea-level at the tide gauge station. This reveals the clear relationship between Ekman transport and sea-level rise.

The simulations associated with this tropical disturbance show that the damaging waves could occur even if the storm is not close to the Senegal coast. The shape of the Senegal coasts makes it vulnerable to strong southerly winds with a long fetch and over a long period. Similarly, [28] show that anti-cyclonic and cyclonic systems are responsible for extreme waves in the British Isles. In Senegal, south of Dakar, developing low pressure associated with African Easterly Waves during the summer can produce southerly winds and thus serve as the source of extreme waves. High pressure can produce northerly winds and would serve as the source of extreme waves in coastal regions to the north of Dakar. These weather events can be predicted with some degree of success, although the prediction of tropical cyclo-genesis is still challenging across the extreme eastern Atlantic Ocean.

Increased Earth observations and predictive weather and wave models provide the best means of understanding events that cause coastal damage along West Africa. Such a system can provide the basis for developing a coastal warning system to protect the millions of people that live in coastal zones in West Africa in large cities, small towns, and fishing villages. The present challenge is to increase in situ coastal data, which has limited real-time buoy data and coastal weather station data. These data are needed to monitor extreme events, evaluate and improve satellite-based products, weather, and wave models.